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I. INTRODUCTION 



Recent experimental developments in heavy ion collision experiments have lead to the 
formation of hadronic matter at finite temperature with high densities. At collisions with 
relativistic heavy ions at facilities like RHIC at Brookhaven and beams at SPS, and eventu- 
ally at LHC, lead to the deconfinement of quarks from the confined state in hadrons forming 
a quark gluon plasma. The phase transition is anticipated at temperature of T c pa 150MeV 
- 250 MeV. The quark gluon plasma above the critical temperature (T c ) and the hot and 
dense hadronic matter below T c possess special problems to investigate their properties at 
finite temperature and high densities. We have to consider interaction of hadrons, the de- 
cay of mesons and neutrons and finally the behavior of quark - quark interactions and the 
screening effect in the quark gluon plasma (QGP). Such a study requires the development 
of a formalism of quantum field theory to deal with problems at finite temperature and high 
densities. 

It should be stressed that, developments to deal with hot and dense nuclear matter and 
hot QGP would provide valuable tools to study early universe as well as the interior of 
stars. Of course there is a range of problems that would require time dependent treatment 
since the system may be in a non - equilibrium state that is changing rapidly with time. In 
general, therefore, it appears that both early Universe and relativistic heavy ion reactions 
would require a formalism that depends on time and temperature simultaneously. 

We wish to explore methods to calculate the reaction rates and decay rates in a medium at 
temperature T using the real time finite temperature field theory - Thermo Field Dynamics 
(TFD) [JTJ] . This theory was developed so that the time degree of freedom is not lost ||. 
It was realized that, an extention of the usual field theory at zero temperature to finite 
temperatures would require doubling of the Hilbert space. At the same time the operators 
have to be doubled such that the new set of operators, designated tilde operators, act on the 
second Hilbert space, tilde space. Later on it was interpreted that, the second Hilbert space 
acts like a heat bath that ensures the dynamical system to stay at constant temperature 
T. It is established that with TFD the Wick's theorem can be applied and other features 
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like Ward - Takahshsi relations and Nambu - Goldstone theorem on spontaneously breaking 
symmetry can be established. Furthermore as in the case of the usual field theory at zero 
temperature, Feynman rules for scattering amplitude and decay rates can be established. 
Therefore the combination of all these feature in TFD allows us to carry out calculations as 
in the case of T = field theory. Some of the essential aspects such as a definition of the 
vacuum and a definition of the creation and annihilation operators at finite temperature are 
given in the Appendix. It should be mentioned that for systems in equilibrium there are two 
other methods available. The first is the imaginary time method due to Matsubara || that 
has been extended to quantum field theory by Jackiw and Dolan Q]. The second method 
is the closed time path method due to Schwinger and Keldysh ||. Both these methods 
are essentually based on Green's function approach. Recently Niegawa has developed a 
technique to use the closed time method to calculate the scattering amplitude and decay 
rates. We'll compare this method with that of TFD later on. 

It should be stressed that what has been done is to write down the decay amplitudes by 
using the method of Cutkosky to find the imaginary part of an amplitude. These rules are 
very helpful in our investigation. 

Recently Jeon and Ellis [6] have developed a multiple scattering expansion of the self 
energy at finite temperature in the imaginary time method. They have calculated the re- 
sponse function in this approach. Brandt et al [6] have calculated the retarded thermal 
Green's function and forward scattering amplitude at two loops. Using time-dependent two- 
and three-point functions in scalar theory they have calculated i) the forward scattering 
amplitude of two on-shell thermal particles and ii) forward scattering amplitude of a single 
on-shell thermal particle. The calculations include one-loop self energy and vertex correc- 
tions that are retarded but are at zero temperature. The present work attempts to calculate 
scattering and decay amplitudes at finite temperatures. 

The paper is organised as follows. In sect. II we briefly present basic TFD formulas, 
which will be nessesary in further calculations. In sect. Ill we describe our basic formalism. 
In sections IV, V and VI we consider decay and scattering procecces respectively. In sect 
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VII we consider loop contributions. The summary and conclusions are given in sect VIII. 



II. THERMAL FIELDS AND S - MATRIX IN TFD. 

The Thermo Field Dynamics is a real time operator formalism of quantum field theory 
at finite temperature. Any physical state can be constructed from a temperture dependant 
vacuum \Q((3)) which is a pure state. The main feature of TFD is that, the thermal average 
of any operator A is equal to its temperature dependent vacuum expectation value with the: 
vacuum, |0(/3)), being obtained from the usual vacuum by a Bogoliubov transformation. 
Therefore, we have 

< A >= (0(/3)\A\0(p)) (2.1) 

where (5 = 1/ksT with ks being the Boltzmann constant. Doubling of degrees of freedom of 
all fields is required. This is achieved through a "tilde" operation: to each zero temperature 
field 4>(x) a doublet of fields (4>(x), 4>{x)) is attached, the dynamics of which is controlled by 
a thermal Lagrangian 

L — L(<f>) — L(<f)) — L(<f>) — L*(<j>) (2.2) 

The temperature dependent vacuum state |0(/3)) is annihilated by the temperature depen- 
dent physical annihilation operators 

0*00)10(0)) = c k {(3)W)) = dMW)) = 

(2.3) 

5*09)10(0)) = c*(0)|O(0)) = d fc (0)|O(0)) = 
which are obtained through Bogoluibov transformation from usual annihilation, a k , and 
creation operators, a\. For a Bose system 

at = cosh.6 k cik{(3) + sinh# fe at(/3), at = cosh^ fe at(0) + sinh 8kCik(/3), 

(2.4) 

a\ = smh9 k a k (/3) + cosh. 9 k a\((3), a k = sinh 9 k a\{(3) + cosh.9 k a k ((3), 

where 

sinh 2 9 k = n B (k) = l/(e^l^°l — 1), cosh fc = e^W 2 sinh fc . (2.5) 
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with ko being the energy associated with the four - vector, k. The creation and annihilation 
operators satisfy ordinary commutation relations 

[a k , aj]_ = [a k , aj]_ = [a k (P), aJ(/3)]_ = [a k {/3), aj(^)l- = 5(k-p) (2.6) 

and all other commutators vanish. 

For fermions, with the creation, c| a and cd, a , and annihilation, c P;S and d p >a operators 
the Bogoliubov transformations lead to the relations 

c PjS = cos 9 +p Cp tS (/3) + ism9 +p c p s (/3), d PtS = cos 6- p d P)S (f3) + i sin d- p d\ >a (/3), 
cj ia = i sin 6 +p c p>s (l3) + cos6> +p cJ a (/3), Jj s = z sin 6_ p d P;S (l3) + cos6L p dJ a (/3), 
4,s = cos^+pCp,^/ 5 ) - «sin6» + pCp iS (/3), c/ ps = cos6L p dJ a (/3) - i sin 0_ p dp )8 (/3), 
c p , s = -i sin +p cj, iS (/3) + cos 9 +p c PtS ((3), d p , s = -i sin 6_ p d) p ^) + cos Q- P d p , s ((3), 
where 

sin 2 ±p = l/(l + e^|P0=FMl), cos0 ±p = e^o±^)/2 sin ^ ±p 
sin 2 +p = n F (p), sin 2 6L p = ^(p) 
Here po is the energy 



[2.7) 



(2.8) 



associated with the four - vector, p and \x is the chemical potential. 
The anticommutation relations for creation and annihilation operators (c- for particles, and 
for antiparticles) are similar to those at zero temperature: 



[ £ W C J,J + = K's' d l,s\+ = [ C pvC#)4,s(/ 3 )]+ = [ d p's'(P) d t,s(P))+ = S(p-J/)5 S 



(2.9) 



All other anti - commutators are zero. In order to define the reaction rates and decay widths, 
we follow the prescription due to Feynman 0, which is best stated by Dyson ||, and this 
involves writing down an operator: 
00 (—i) n r+oo 

H F (x ) = J2 ^-h / dx x dx 2 ■ ■ ■ dx n T[H e (x )H I (x 1 )H I (x 2 ) . . . H I {x n )\, (2.10) 
n=o nl J -°° 

Then choosing H e (xo) = 1 this reduces to an expression that is evaluated between an initial 
and a final state. Therefore we get an expression that looks like and is equal to an S matrix 
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00 00 (—i) n r 

S =Y,S n = J2 / dx i dx 2 ■ ■ ■ dx n T[Hi{xx)H I {x 2 ) . . . Hj(x n )}, (2.11) 

n=0 n=0 n - 

For the total transition in TFD we have 
00 00 (—i) n r 

S=J2S n = J2 / <bidx2 ■ ■ ■ dxnTlHjix^H^) . . . ffj(i n )], (2.12) 

n=0 n=0 n - 

where Hi(x) = Hi(x) — Hi(x) 0. Particularly, at the tree level 

Sttl-i J d i x(H I (x) - H I (x)) (2.13) 

Now in these expressions Hi{x) and Hj(x) may be introduced explicitly and then Wick's 
theorem is used as in the case of T = field theory to obtain all the contributions in 
perturbation theory. The evaluation proceeds by using a set of Feynman rules for evaluation 
of a reaction rate and decay widths. It is clear that for operators appearing in Hj(x) and 
Hj(x), the transformation of boson and Fermion operators given by Eq.s (|2.4j) and (|2.7|) 
respectively has to be used to bring in the temperature dependent factors. The one particle 
Green's functions are given explicitly in the Appendix. 

III. CALCULATION OF REACTION RATES 

Let us consider the process 

Pl +p 2 + ■ ■ ■ +p n -»• P'l +P2 + ■ ■ ■ +P'm (3-1) 



At zero temperatute, T = 0, the amplitude of this process can be calculated by usual 
Feynman rules: 

oo 

(f\S\i) = j:(f\S^\i) (3.2) 

n=0 

where |z) = a^a^ . . . a|JO), |/) = a pl a p , . . . a} p , |0) and |0) is an ordinary vacuum state: 
a p \0) = 0. In the present article we propose that the amplitude and, hence the cross section, 
of process (|3.1|) can be directly calculated perturbatively using Eq.s ( |2.12| ) at T ^ also. 
Namely, we propose that 
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1) The transition amplitudes are determined from Eq. (|2.12|) where for the process ( |3 . 1| ) 
we have \i) = a^a^) . . . 4 n (/?)|0(/3)) and |/) = 0^(^)0^/3) . . . 4^(/3)|0(/3)) . 

2) The phase space factors, relating an amplitude to the rate are the same as in zero 
temperature case. For example, the differential cross section for the process pi + P2 —>■ 
p\+p' 2 + ...+p' m is 

d " = < 2 ^ 4 « +14 + ... +14. - - n< 2 '"'> ft pf^fil 2 < 3 - 3 > 

where _E P = \jra 2 + (p) 2 and t> re j is the relative velocity of the two initial particles with 
momenta p\ and p 2 - The amplitude is related to S matrix in manner similar to the procedure 
carried out at zero temperature 

(f\S\i) = i(2^5\P f - Pi )M & nt^] 1/2 n(^] 1/2 - (3-4) 

ext Vtj ext VUJ 

Here Pi and Pj are the total four momenta in the initial and final states, and the products 
extend over all external fermions and bosons, E and us being the energies of the individual 
external fermions and bosons respectively, and V is the volume. We propose that, in TFD 
this definition still holds. 



IV. DECAY OF PARTICLES 



We know that at T = the decay rates T can be calculated either by studying disconti- 



nuities of the self energy of a decaying particle (Cutkosky rules [[yj]) or by direct calculation 
of the amplitude perturbatively (Feynman rules). Obviously these two methods must lead 
to the same result. 



For finite temperature the Cutkosky rules were generalized by Kobes and Semenoff [JT2 
In the next section we shall calculate the decay of a particle at finite temperature 

a) by using generalized Cutkosky rules (GCR); and also 

b) directly as a square of module of elements of S matrix, i.e. by using the above method. 
It will be established that, these two methods give the same result for T. We shall present 
these to establish their equivalence. 



A. a(k) -> tt°(Asi) + n (k 2 ) 



Firstly, we calculate the decay rate of sigma meson into two pions (we'll omit isospin 
indices for simplicity) in real time formalism, using GCR. The interaction Lagrangian is 
-^int = ^ a ( x ) Tc2 ( x )- The decay rate of a boson with mass M and four momentum k = (a>, k) 



with u) = Y M 2 + (k) , is related to the self energy [III] by 



T GC r(w) = --ImE(u) = - 1} /mSn(fc) (4.1) 

w w(eP^ + 1) 

where 

iSu(Ar) = A 2 / ^iAS^iAkCp- *) (4.2) 

The GCR fl^], as illustrated in Fig. 1 give: 

2ImE u (k) = X 2 J -^^[iA + (p)iA-(p-k)+iA-(p)iA + (p-k)] (4.3) 

where iA (p) = 27r[@(±po) + nB{p)}5(p 2 — m 2 ). The two terms in the integrand are related 
by 

iA ± (p)=ie ± PP°A*(p) (4.4) 

from flO]) - (|PD we get: 

n / \ ieP u -l)(e-^ + 1)A 2 /• dp 4 , 2jf<r 2 3 ,_„ n2 2 , 
r G cflM = ^ /V- — / -f—A2Ti) 2 5(p 2 - m 2 )8((p - k) 2 - m 2 )x 



2u(e^ + l) J (2tt) 4V yy ^ (4.5) 

x[©(Po) + tib(po)][©(-Po + + riij(po - 

We shall proceed and evaluate r^oR explicitly. Note, however that, at finite temperature 
Lorentz invariance is lost, and hence, the decay rate will no longer be Lorentz invariant, but 
dependent on the reference frame. In the following examples we will choose the rest frame 
of the decaying particle, as the reference frame, that is we set k = 0. First, we consider the 
product of the mass shell S functions. With k = (M, 0, 0, 0) the compatible zeros are easily 
found to be p = u p = M/2. Hence the delta functions reduce to : 
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5(p 2 - m 2 )6((p - k) 2 - m 2 ) = j^j-JiPo ~ uj p )6(^p ~ M/2) (4.6) 

and fix the momentum dependence of the integrand completely. 

The thermal factors in (|4.5| ) can be reduced on mass shell: po = M/2, u = M 
( e f3u _ i)( e -(3u + i) 

; n) 46(po) + n B ( Po )][Q(-p + u>) + n B ( Po - u)] = 

{eP u + 1) (4.7) 

= (1 _ e -fa)[i + n B (po)][l + n B ( Po - W )] = (1 - e~P M ){\ + n B (M/2)) 2 



now inserting (4.6), ( 4/7) into (fOj) we get: 
T GCR (T^0) (l + n B (M/2)) 2 



'l + n B {M/2)f -n 2 B {M/2) (4.8) 



(4.9) 



T(T = 0) l + n B (M) 

Here we used the simple relations: 

e -(3M = e -(3M/2 e -(3M/2 

e~P u = n B (u)/{l + n B (u)) 

which will be useful for further calculations. 

Now we calculate the rate T(a — > im) directly. Because of the doubling of degrees of 
freedom the whole interaction Lagrangian will be Lj n ^ = \a(x)7r 2 (x) — \a(x)7r 2 (x). The 
initial and final states are \i) = a k (/3)\0((3)) and |/) = b\. ((3)b\. (/3)\0((3)) where operators 
a((3) and b{(3) stand for a and 7r fields respectively. At the tree level (??) the matrix element 
of S matrix is given by 

(f\S\i) =i\J dx(0((3)\b k MK(f3)lv(^ 2 (x) - a(x)n 2 (x)}a{([3)\0({3)) (4.10) 

where the boson fields a(x) and tc(x) are defined as in eq. ( |A22j ). Using Bogoluibov trans- 
formations ( |2.4j ) and commutation relations ( |2.6j ) it is easy to show that: 

(0(P)\a(x)al(P)\0((3)) = e~ ikx cosh 9 k , 

(0(J3)\*(x)4(P)W)) = e- ikx smh9 k , 
W)\b kl {f3)b k2 {f3)it 2 {x)W)) = jfa + ^cosh^ cosh^ 2 
(0(P)\b kl (P)bk 2 m 2 ^)W)) = e^ 1 + k2 >smhe kl sinh^ 2 . 
/From Eqs (|4.10| ), ( [4.11] ) and (|2.12|) we obtain the amplitude: 



(4.11) 



Ai^iT) = A [cosh 6k cosh 9^ cosh^ — sinh 9k sinh 9^ sinh^ 2 ]. (4-12) 



The decay rate, can be calculated directly from (|4.12j ) without introducing any additional 
thermal factor: 

1 r dkxdk 2 {2Tt)%k-kx-k 2 )\M§{T)\ 2 



2ujJ (2cj 1 )(2tu 2 )(27r) 3 (27r) 3 
A 2 f dkxdk 2 S 4 (k - h - k 2 )W B (u, wi, u 2 ) _ X 2 I B (T) 



(4.13) 



32cj7t 2 J u)\U)2 32um 2 

where 

Wb(w, Mi, UJ2) = [cosh Ok cosh 6'fc 1 cosh 9k 2 — sinh ^ sinh 9k x sinh #/c 2 ] 2 , (4-14) 



= y (fcj) + m 2 and = Vv^) + M 2 . Now using the properties of the hyperbolic func- 



tions and Eqs ( |2.5|) , and ( [4.9| ) this may be further simplified as follows: 

5(uo -ui- w 2 )W B (w, u h u 2 ) = n^n^el 3 ^ + Ul + ^ 2 )/ 2 - l\ 2 5(u -u x - u 2 ) = 

= n x n 2 n^e^ - \){eP^ + U2 - - wi - uj 2 ) = 

1 + 1%, (l + m)(l + n 2 ) , ( ^ 4 ' 15 ' ) 

= rii^n^l 1}{ l}d(uj -u x - ui 2 ) = 

fix n 2 

= (1 + nx + n 2 )5(uj — ojx — uj 2 ) 

where for simplicity we denote = n B (ki), (i = 1,2), and = n B {k). 

The integral I B {T) in Eq. ( 4.13| ) can be explicitly calculated in the rest frame of the 



decaying particle: to = M, k = 0, uoi = y (A;*) + m 2 = J (q) 2 + m 2 = uo q : 
dk 1 dk 2 5 4: (k — k x — k 2 )W B (to, Ux, u) 2 ) 



UJlUJ 2 

^ d qq H^ q -M)W B {M^ q ) = 8n / m _ {m/M y WB{MM/2M/2) 

UJ 2 

(4.16) 



using ( |4.15| ) and ( fl.lb] ) in ( |4.13| ) we finally obtain: 
T(T ^ 0) 



W B (M, M/2, M/2) = (1 + 2n B {M/2)) (4.17) 



r(r = 0) 

this formulae is the same as obtained by the GCR method eq 
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B. H(k)^e-(h) + e+(k 2 ) 



As a second example we consider the decay of a boson into two fermions, namely, the 
decay of Higgs boson into electron - positron pair. The decay rate of this process has been 



calculated in detail using GCR by Keil 14]. The final tree level result is : 



T GC r(T I 0) _ (1 - n F (M/2))(l - n F (M/2)) 
T(T = 0) l + n B (M) 1 ' ] 

which is obtained in the rest frame of H boson with a mass M. 

Now we evaluate the same rate by direct calculation of the amplitude. The interaction 
Lagrangian is 

Ant = L int - Ant' L int = -igH(x)ij){x)if>(x), L int = +igH(x)ip(x)tp(x). (4.19) 

The matrix element of S matrix from an initial state \i) = a k (/3)\0((3)) to the final state 
\f) = 409)4,00)1000)) at the tree level is 



(f\S\i) = i(-ig) I dx 4 (i\[H(x)^(x)tlj(x) + H(x)ip(x)ip(x)]\f) (4.20) 
where we omit the spin indices for simplicity. As in previous example, Eq. ( |4.11|) , we have 

(/|#(a;)|i) = e~ lkx cosh^, {f\H(x)\t} = e~ ikx smh9 k . (4.21) 
For the fermions, using ( |A23| ) in ([4.20 ) we have 



(f$(x)iP{x)\i) = J dpdfiN p N' p (0(P)\c k MdkMHHpy pX + d^e-^x 
x [c p/ u(p')e- i P' x + d^vip'yV^W)) 

and 

(f$(x)j;(x)\i) = J d^iVpJV^O^Icfc^^dfe^^fi^e^ + ^^e^X 
x [c pl u{p')e l P x + 4^(p')e~^ /x ]|0( / 9)) 

The only "surviving" term in ( 4.22Q is that one which includes 
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(4.22) 



(4.23) 



(0(P)\c kl (P)d k2 (P)4dlW)) = (0(i3)\c kl (l3)d k M[^se +p cl(i3) -isme+pCpiP)] = 

= [cos^4,(/3) -zsin^^(/3)]|0(/3)> = (0 (/3) | c fel (/3) d fe2 (/?) c+ (/3) 4, (/3) 1 0(/3) > x (4.24) 
x cos 9 +p cos 9-pi cos9 +kl cos 9- k2 8(k\ — p)S(k 2 — p'). 
Similarly, the nonzero contribution to ( 4.23|) involves 

(0(f3)\c k MdkMd P 5 P >\0((3)) = (0((3)\c k MdkM[cos9_ p d p ((3) - isiYi6„Mf3)}* 



x [cos 9 +p ,c p , {(3) - 1 sin +p ,cL ^)]W)) 



(4.25) 



= - (0(/3) \c kl (f3)d k2 (f3)cl, (P)dl(P) |0(/3)> sin 9_ p sin 6 +p , 
= - sin6» +fcl sin6»_ fc2 (5(A;i - p')5(k 2 -p). 

The decay of the process is given by a general formula of the zero temperature field 
theory: 

1 r dhdk 2 {2^\k ~h- k 2 )(2m) 2 \M^(T)\ 2 
[UJ) 2uJ (2cui)(2^ 2 )(27r) 3 (27r) 3 

where the transition amplitude is given by Eqs (|3.4|) , and for the present case it has the form 



A4ft(T) = (—ig)[cos9 kl cos9- k2 cosh 9 UJ u(ki)v (k 2 ) — sin 6^ sin9- k2 sinh 9 u v {k 2 )u{k\)]. (4.27) 



Then Eq. ( f4.26| ) can be fiurther simplified by using 



|.Mg(T)| 2 = g 2 [cos9 kl cos 9- k2 cosh 9 W — sin 6^ sin 9- k2 sinh 9j\ 2 Tr 



{jk 1 + m)(jk 2 + rn) 



spins 



2m 



2m 



(4.28) 



and 



5{oj — lji — u 2 ) [cos 9 kl cos 9-k 2 cosh 9 U — sin 9 kl sin 9_ k2 sinh 9 U ] 2 = 

= 5(uj - wi - c^ 2 ) cos 2 9 kl cos 2 #_ fc2 cosh 2 W (1 - e~P( w + ^ + ^a)/ 2 ) 2 = 

= 5{uj — u)\ — co> 2 )(l — np{k\){\ — n^(A;2)(l — e^^ u ) = 5(u — oj\ — uj 2 )Wf(oj, u>i, uj 2 ). 

(4.29) 



Now, inserting (|4.271 ) - ([4.29|) into (|4 . 2 6|) and performing integration in the rest frame of the 
decaying boson we obtain 
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r(T ^ o) 



W F (M, M/2, M/2) = (1 - n F (M/2))(l - n F (M/2))(l - e~ (3M ) (4.30) 



T(T = 0) 

which is exactly the same as Eq. ( [4.18| ) given by GCR, since (1 — e~@ UJ ) = (cosh 2 ^) -1 



[l + n B (uj)) 



-i 



C. $(k) ^ <P(h) + <P(k 2 ) + <j>(k 3 ) 



The next example is the decay rate of a boson -$(a;) into three bosons - <f)(x) with the 
interaction Lagrangian = — — &(x)(J) 3 (x), and L int = — L int . . The decay rate of 
the process (both in imaginary time formalism and in TFD) using GCR has been obtained 
by Fujimoto et al. |13|] to be 

ImTj(uj) f dk\dk 2 5{uj — uo\ — uo 2 ~ uj 3 ) 

GCR\ ~~ 
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«0 = — = ^ o\bo . . " >< 



(27r) 6 8^iu;2W3 (4.31) 



x [(1 + m)(l +n 2 )(l + n 3 ) - nin 2 n 3 ] 



where n; = n F (fcj), Ui = y (fcj) + m 2 , (i = 1, 2, 3) and co> = y(A?) + M 2 . On the other hand 



the decay rate is related to the transition amplitude as: 

1 r dhdhdk^n^djk-h-h-h^M^T)] 2 
{UJ) 2uJ (2 Wl )(27r) 3 (2o; 2 )(27r) 3 (2o;3)(27r) 3 

where the amplitude Ai^(T) is defined by Eq. (|3.4j ) with the elements of S - matrix: 

(f\S\i) dx{m\K^)KWkM[^)^^) ~ Hx)4> 3 (x)]a{((3)\0((3)} (4.33) 



As in the previous example, we get 
5(u — oj\ — u 2 — uJ S )\Mj^(T)\ 2 — 

= [cosh 9% cosh 9^ cosh 9k 2 cosh 9k 3 — sinh 9k sinh 9^ sinh 9^ 2 sinh 9k 3 ] 2 S(uj — uj\ — uj 2 — uj 3 ) = 
= ni n 2 n^[l - eP( u + ^1 + ^2 + uj 3 )/2^ 5 ^ _ ^ _ ^ _ ^ = 

= nin 2 n 3 n £J [l - e^][l - + ^2 + ]«$(<,, - u; 2 - cj 3 ) 

1 + (l + m)(l + n 2 )(l+n 3 ) 

= n x n 2 n z n UJ 1 1 \d(uj - Ui - uj 2 - lo 3 ) = 

riu nin 2 n 3 

= [(1 + ni)(l + n 2 )(l + n 3 ) - nin 2 n 3 ]5(u - u) X - u 2 - u 3 ) 

(4.34) 
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where n w = riB(k), n« = Unifa). Then, by using ( f4.34|) in ( |4.32| ) the result given in Eq. 



( |4.31| ) is reproduced. Again the direct calculations and those from GCR agree. 



V. THE DETAILED BALANCE. 

The last example of the previous section can be easily generalized for the decay of a 
boson into m particles: <3>(fc) — > 4>{k,i) + 4>{k2) + . . . + 4>{k m ) with the decay rate 

_ }_ f dhdfo • • • dk m {2iiY5\k - h - k 2 - . . . - k m )\M^{T)\ 2 
^' 2uJ (2u 1 )(2u< 2 )...(2u m )(27i) 3m 



where 



(5.2) 



6(U) -Ui - . . . - U m )\Mft(T)\ 2 = 5(UJ - U)t - . . . - LUm) 

[cosh 9k cosh 9^ . . . cosh 9^ m — sinh 9^ sinh 9^ . . . sinh 9k m } 2 = 
= n x . . . n m n w [l - eP u \ [1 - + ■■■ + ^m)] = 

= [(1 + n x )(l + n 2 ) . . . (1 + n m ) - n x n 2 . . . n m }5(u - u x - . . . - u m ). 

It is clear from Eq. (|5.1|) and ( |5.2|) that, the total decay rate T(u) is the difference of 
the forward rate of the boson decay - T^{oj) and the rate of the inverse process Yi(uj): 
T{uj) = r d (uj) — Ti(u), which are related by the principle of detailed balance: 

= (l + ni)(l + 7i 2 )...(l + 7i m ) = + U2 + . . . + Um) = e(3u 

VAu) n x n 2 ...n m 



This relation was first shown by Wei don, [T^| by analysing the discontiniuties of the self 
energy function in the Matsubara formalism. Thus, we conclude that, our approach of direct 
calculation of rates through elements of the S matrix leads to the correct relation for the 



detailed balance principle at finite temperature. Note that, the Eqs. ( |5.1| ) - ( |5.3|) reveal the 
importance of doubling of degrees of freedom at nonzero temperature: If we had neglected 
the term then the term proportional to "sinh^" in Eq. fl5.2|) would have vanished 
giving a wrong relation. Although the "tilde" particles are introduced as an essential part 
of the finite temperature formalism as fictitious particles, they play an important role in 
direct and inverse processes taking place at finite temperature. It is important to remark 
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that there are no transitions, and hence no matrix element, between the tilde and non - tilde 
particles. 

However, following Keil [Q, we underline that, T(u) and not the partial rates Td(uj) and 
Ti(u) represent the physically measurable decay rate! The reason is that, it is the full decay 
width T(u) = T d (uj) — T^Ld) that is connected to the pole of the Green's function |13|,?1 

l + n B (k ) n B (k ) 



An(fc) 



k 2 — ml — £ + ie k 2 — m 2 , — E* — ie 
l + n B (k ) n B (k ) 



(5.4) 



(fc + ^/2) 2 - uj 2 + ie (k - tT/2) 2 - oo 2 - ie 



where uo 2 = k 2 + ReY* - T 2 /4, V = -ImE/k , k = yj(kf + m 2 . 



VI. SCATTERING CROSS SECTION OFl + 2^1' + 2'. 

As it was stated in Sect. Ill in our approach, in accordance with the concept of TFD, 
the relation between the rate (cross section) and the transition amplitude is the same as it is 
in the zero temperature quantum field theory. For example the cross section of the process 
a(h) + b(k 2 ) -> a'(ki') + b'(k 2 ') is [0 

da _ |^ fi (T)| 2 (n ; (2mQ)|fcT| 2 ^K + ^) 
dVl' Q4Tc 2 v Te \uJiUJ2Uj'iUJ2 d\ki'\ 

in the usual notation. In particular, the cross section of the elastic scattering of two bosons 

in their center of mass (CoM) frame [k\ = —k 2 ) is 

da_ \M R (T)\ 2 
W 647T 2 (cu 1 +o; 2 ) 2 

These equations imply that 

(da/<KV)\ T # _ \M^T^0)\ 2 
(da/dn')\ T=0 \Mz(T = 0)\ 2 

where the amplitudes may be determined through Eq.s (|3.2j )- (|3.4| ). To see the consequences 
of this relation we will consider some examples. 

A. Boson - Boson scattering. 



(6.2) 



= W(T) (6.3) 
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Let us assume that both a and b particles are bosons with the interaction Lagrangian 
-^int = ^a( x )^b{ x ) an d L int = -^int — ^int ■ At the tree level, using the results given in 
Eq.s Q - O and Eq. (pi]) in Q we get 



l-Wfi(T = 0)| 2 ( 64 ) 
= m^^^^i + ^ + wi + <4)/2) _ l]2 j 



where 



(6.5) 



C(T) = cosh 6^ cosh ^ cosh^' cosh# fc2 ', 
S'(T) = sinh 6^ sinh 9 k2 sinh 6* fcl ' sinh 9 k2 ', 
and = ne(fcj), ^ = Due to energy conservation we have 

[ e /3(wi + uj 2 + uj[ + J 2 )/2 _ ^2 = ( e /3(^i + w 2 ) _ !)( e /3K + Wa) _ i). (6.6) 

Now we may express the exponents via n^, n[ 

= I±*, = (6.7) 



to get 



W flB (T) = (1 + m + n 2 )(l + n[ + n' 2 ) = 

= [(1 +ni)(l + n 2 ) -nin 2 ][(l + 0(l +n' 2 ) - n> 2 ] 



(6.8) 



B. Fermion - Fermion scattering. 

Now let us assume that both a and b particles are fermions with the interaction la- 
grangian Lj n |. = gip a (x)T ol ip a (x)ipb{x)r a ipf } (x) where T a is an appropriate combination of 
Dirac matrices. Similarly as in the previous case we get : 

W FF {T) = !^ fi ^ - = [cos 9 kl cos 9 k , 2 cos 9 kl > cos 9 k2 > - sin 9 kl sin 9 k , 2 sin 9 kl > sin 9 k2 >} 2 = 
\Mfi{l = U)| 

= n F (A; 1 )n F (fc 2 )n F (fc 1 , )n F (A; 2 , )[e^(^ + w 2 + £«4 + w 2 ~ Ml ~ A*2 - m'i - /i' 2 )/2 _ 1 j 2> 

(6.9) 

At equilibrum /ii + /i 2 = // x + /i' 2 , and hence, by using 
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[ e P(vi - Hi + u 2 - _ i]2 = [ e /3(^i - /ii + w 2 - ^2) _ 1] [e/3(<*>i - A*i + ^'2 ~ /4) - l] 
1 1 1 „ 1 1 1 , 



n F (ki)n F (k 2 ) n F (k 1 ) n F (k 2 ) n F (ki)n F (k 2 ) n F (kx) n F (k 2 ')" 



(6.10) 



we get 



W FF (T) = = [(1 - ni )(l - n 2 ) - mn 2 ][(l - ni)(l - n 2 ) - n^]. (6.11) 

where = n F (ki), and raj = n F (ki). 
C. Fermion - Boson scattering. 

Let us assume that particle a is a fermion and particle b is a boson with the interaction 
lagrangian given as = —igip(x)Tip(x)(f)(x) ( for example 7 = 75 for pion - nucleon 
scattering.) In this case we get at tree level 

W FB (T) = |"^ fi ^ ~ °2\ 2 = [cos 9 kl cosh 9 k2 cos 9 kl > cosh 9 k2 > + sin 9 kl sinh 9 k2 sin 9 kl > sinh fca /] 

(6.12) 

Then using similar manipulations, as in previous examples we get 

W FB (T) = l^j^yp = [1 " Mfc) + n B {k 2 )][l - n F (h') + n fl (k')]. (6.13) 



The expressions (|6.8|) , (|6.11|) and ( |6 . 13|) for the relation of in medium cross sections of 



1 + 2 — > V + 2' scattering can be interpreted physically in the spirit of Weldon |15|]. For 
example W F b(T) can be rewritten as 

W FB (T) = (1 - n F (h))(l + n B (k 2 ))(l - n F (A; 1 / ))(l + n B (A: 2 '))+ 

+r^ F (A; 1 )r^ i J(A;2)n F (A; 1 , )n B (A;2 , ) + nj^OMfeX 1 - n F (fci'))(l + n B (k 2 ))+ (6.14) 

+np(A;i')nB(A; 2 ')(l - n F (h))(l + n B {k 2 )). 

It is to be emphasized that the in-medium cross sections at finite temperature are not directly 
measureable in any many-body process. However such cross-sections are an important input 
to treat a gas at low density through a kinetic theory approach such as a Boltzmann equation. 
The finite temperature cross-sections will be a part of the collision term that would lead 
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to a distribution function by a solution of the Boltzmann equation. As indicated in the 
introduction these types of results are a necessary input to a better understanding of both 
the early Universe and the hadronic matter at finite temperature and density. 



VII. LOOP CORRECTIONS 

All the above results have been obtained at the tree level. Now we wish to show how the 
method can be handled beyond the tree level. For simplicity we choose the (^(x) interaction 
•-^int = ~ : 4 ( x ) : • The S matrix in Eq. ( |2.12| ) up to order A 2 is given by: 

S^l-tJ dxH^x) + J dxdyTlH^H^y)} = S {1) + S {2) . (7.1) 

For the scattering process from the initial state 

|<> = = alW)al(P)W)) (7.2) 

to a final state 

\f) = mMh')) = g^(/3)4 2 ,(/3)|0(/3)) (7.3) 

the matrix element (f\S^\i) and hence the amplitude A4^(T) = M.^(T) + .Mg^T), is 
defined as, 

(f\SW + S^\i) = 5 if + 5\P - P') fi fi - , (7.4) 

(An) 2 J ujiuj2^i^2 

with P = h + k 2 and P' = ki + k 2 '. Mf{T) and M^iT) can be simply calculated by 
using the procedure used in sec. VI. For instance, 

M ( l\T) = -zX[C(T)-S(T)} (7.5) 



where C(T) and S(T) are defined by Eqs 

To evaluate Ai^ (T), the generalized Wick theorem [ To | is used and 2x2 Green's function 
A a b(q) (see Appendix) is introduced. In particular: 
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iA° 22 {x - 



y) = W)\T[$( X )$(y)]\0((3)) = i J ^e~^ x ~ ^A° 22 (g) (7.6) 



d A q 

In terms of these Green's functions, the amplitude is proportional to A 2 and is 



zA? 2 (x-y) = (O03)|T[^(y)]|O03)> = </ ~ v) A° 2 (g). 



M$(T) = A 2 {C(T) J ^- 4 A u (q)A° u (P - g) + S(T) J ^A° 22 (q)A° 22 (P - q) 
dq 4 



(2tt) 4 UV1 ' iiV v y 7 (27 

(2 ^ )4 [C s (/ci, fc 2 ; An', fc 2 ')A° 2 (g)A? 2 (P - g) + C s (h\ k 2 '; h, k 2 )A° 21 (q)A° 21 (P - g)]}. 

(7.7) 

where C s (k±, k 2 ; k±, k 2 ) = cosh 9 kl cosh^ sinh# fcl ' smh9 k2 >. The terms in the integrand are 
related by 



AM?) = -[A^g)]* = ,_ 2 - - 2^ B (?)5(? 2 - m 2 ), 

y " t i ffc (7.o) 

A? 2 (g) = A^g) 
and hence Eq. (|7.7f ) simplifies to 

Mj?'(T) = A 2 {C(T)P(P,T) + 5(T)F*(P,T) + 

+87r 2 e^/ 2 / 3 ( J P)[C s (fc 1 , fc 2 ; fc/, fc 2 ') + C s (ki, k 2 '; h, k 2 ))}, (7.9) 
F(P ' T) = / (2^) A?l(g)A?l(P " g) = /l(P) ~ 47U/2(P) ~ 4?r2/3(P) 



with 



dg 4 1 



(2tt) 4 (g 2 - m 2 + ze)((P - g) 2 - m 2 + is) 



The integrals I\(P) , I 2 (P) , I^(P) can be evaluated in the CoM frame of the two particles , 
where P 2 = (k® + k 2 ) 2 = E 2 - total energy of the system. 
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The integral l\ [P) is logarithmically divergent and can be treated using a renormalization 
procedure and this finite contribution is called the vacuum fluctuation correction. This term 
exists even at zero temperature, and can be calculated by the usual method of dimensional 



regularization | 17fl . As a result h(P) may be written as: 



iN F i r , , m , vfl 



h(P) = -r^ + tt-ME) -ln^- 



167T 2 167T 2 fl 2 ' 167T 



6o(^) = 2 + v A^ln|i^^| (7.H) 



D = 1 - 4m 2 /£, 

where the first term involving N e = 1/e — 70 + ln47r is infinite. To eliminate the latter we 



introduce the following counterterm in the Lagrangian [17 



\2iu 

Wr = "ig^jfj : * 4 (*) : (7-12) 

which clearly does not depend on temperature. On the other hand, in accordance with the 
tilde oparation rule, this automatically leads us to introduce its tilde partner also: 

A 2 iV 

4ontr = -T^fi : ^ : ( 7 - 13 ) 
which serves to compensate the divergence, coming from S(T) * F*(P, T) term in Eq. ( [7.9|) . 

The ^{P) is convergent and can be evaluated numerically: 

dq f+°° dq n B (q )[5(q + \w q \) - 5(q - \u q \)] 



UP) 



(27r) 4 y-oo 2^((P-g) 2 -m 2 ) 
dqq 2 n B (q) _ 1 



(7.14) 



r/ 2 (£,T) 



4tt 3 J E 2 (l - Au 2 /E 2 )u q 64tt 3 

The integral ^3 including two 5— functions may be evaluated in a similar way as in sect. 

IV (Eq. (|4.6| )) and is given as 

tcp\ [ dq dqn B (P - q)n B (q) n 2 B (E/2)y/D ~ 

HP) = J (2^4B^ 6{q ° ~ ^ K ~ E,2) = • (7 - 15) 

Now , using Eqs. ( |7.5|) , (|7.9|) - (|7.15|) and introducing A^g = iAi^ gives 
A7 fi (T) = ReMftiT) + UmM^T) 

ReM & (T) = A(C(T) - S(T))(1 - XA(E,T)) (7.16) 
ImM^T) = -A 2 v^[(C(T) + S(T))(1 + 2n|) - 4C s (T)e^ 2 n 2 B ]/(16n), 
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where A(E,T)=b (E)-\n^-I 2 (E,T),n B = n B (E/2),C(T) = (l + n B ) 2 and S(T) = n%, 
C S (T) = 2e /3E / 2 n B and A = A/16tt 2 

From ( |7.16|) one may get the relation for the in - medium elastic cross section in 4> 4 (x) 
scalar theory up to one loop approximation: 

(da/<KV)\ T¥0 _ (C(T) - g(T)) 2 (l - XA(E, T)f + 7r 2 A 2 D [(l + 2n|)(C(T) + S(T)) - 8eP E n%] 2 
(da/dM) | T =o ~ [1 - \{b {E) - In ^)] 2 + n 2 D \ 2 

(7.17) 

The appropriate Feynman rules for <f) A scalar theory at finite temperature can be ex- 
plicitly obtained from Fig. 2. They are similar to those for zero temperature. The main 
differences are the following. 

The amplitude of any order n consists of n vertices. There are two kinds of vertices: "dot- 
ted" (ordinary) and "circled" ones. To calculate A^g(T) 

1) Attach factor — iX to each "dotted" vertex. 

2) Attach the factor cosh# p to each external leg (incoming or outgoing particle) with mo- 
mentum p connected to the "dotted" vertex. 

3) Attach the factor sinh^p to each external leg (incoming or outgoing particle) with mo- 
mentum p connected to the "circled " vertex. 

4) Attach iAn(g) to the propagator carrying momentum q, connecting two "dotted" ver- 
tices. 

5) Attach zA 2 2(<?) to the propagator carrying momentum q, connecting two "circled" ver- 
tices. 

6) Attach iAi2(q) to the propagator carrying momentum q, between a "dotted" and a "cir- 
cled" vertices. 

7) Attach iA 2 i(q) to the propagator carrying momentum q, between a "circled" and a "dot- 
ted" vertices. 

Using these rules one may calculate any process in 4 theory at finite temperature to any 
order of perturbative theory. 

Gluon Self-energy and Screening length. 

The gluon self-energy in one looop order is calculated in order to obtain the screening 
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length in a quark-gluon plasma. The calculation is straight forward in Thermo Field Dy- 
namics. Working in Feynman gauge the ghost states do not contribute. Then the photon 
self-energy ir(k, ko] T) in the limit k — > and k — > but at finite temperature reduces to 
a simple form. The screening length is related to real part of the self energy which has the 
expression in one loop order as 

Renf (k - 0, k - 0; T ± 0) = ^-^(N + ±JV»r(2)£(2) = fi*±(N + ^% 2 T 2 

where Nf = number of flavors and a, b refer to the color indices of the group SU(N). 
Here the term proportional to N arises from the triplet and quartet couplings of gluons and 
the term proportional to Nf arises from quark-anti-quark loops. Then the color electric 
mass of the gluon is 

m 2 ee = 1 -(N+ 1 -N f )g^ 
and the screening length at finite temperature is 

L(T) = . £ 

+ \N f gT 

The color magnetic mass of a gluon at finite temperature is zero. Details of this calculation 
are given by Zhao and Khanna [23]. 

VIII. DISCUSSIONS AND SUMMARY 

In the present paper the decay amplitudes and scattering processes are studied in TFD. 
We have proved that, in the TFD formalism, the rate and transition amplitude of a pro- 
cess can be evaluated in a similar way as in the zero temperature case. Note that, both 
imaginary and real parts the amplitudes at finite temperature can be calculated directly. 
The amplitudes are written by establishing a set of Feynman rules using the Lagrangian 
L = L — L. The Feynman diagrams are considered. Use of Bogoluibov transformations to 
define a pure vacuum state leads to Feynman rules for including temperature dependant 
factors for a process occuring in a system at finite temperature. It will be interesting to 
evaluate for example iVW scattering in nuclear matter at finite temperature. Temperature 
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dependence of coupling constants and masses have been explored in perturbative theory of 
TFD [|L8|,?1. At present we are attempting a self consistent calculation of these quantities. 
The role of Ward Takahashi relations at finite T will be crucial to obtaining a self consistent 
solution. 

The present discussion for nuclear matter in equilibrium at finite T can be easily extended 
to QGP at temperature above the deconfinement phase transition. Some of the interesting 
properties like screening length and quark - quark interaction at finite temperature have 
been investigated at finite T |2(J. The results are obtained in a simple and straightforward 
manner. There are numerous other properties of nuclear matter and QGP that may be 
calculated and would be helpful in studies of many particle hadronic systems. 

The method developed in this paper using Thermo Field Dynamics are applicable to 
calculations at tree level as well as loop calculations. This will allow us to calculate processes 
in nuclear matter as well as in quark gluon plasma at finite temperature. 

The present approach simplifies calculation of concrete processes taking place in hot 
and dense matter. Examples are neutrino electron scattering in thermal plasma pi,?] , (3 



- decay, which plays an important role in cosmology, in -medium nucleon - nucleon cross 
sections |23] taking place in ion - ion collisions etc. These processes are being investigated 
with the present approach. 
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APPENDIX A: 



A. The nessessity of doubling. 

In the statistical mechanics the statistical average of a quantity A is given by 

< A >= Z- l {(3) r Ti[Ae-f 3H } (Al) 

where H = H - /iN , Z((3) = Tr[e - ^]. In the zero temperature field theory the average 
of an observable A is determined by it's vaccuum expectation value: 

<A>=(0\A\0). (A2) 

To make the temperature formalism completely parallel to the zero temperature case we 
assume that there exists a " thermal vaccuum state" , such that 

< A >= (0(P)\A\0(J3)) = Z- l J2e~^ En {n\A\n) (A3) 

n 

for any operator A. Since the Hilbert space is complete we can expand |0(/5)) in terms of 
\n) as 

W)) =£|n>(niO(/?)> =Y,fM\n). (A4) 

n n 

^From the last equation we obtain 

W)\A\0(f3)) = f*MUP){n\A\m) (A5) 

n,m 

which will agree with eq. (|A3|) provided 

fn(P)fm(P) = Z^e'^H^. (A6) 

However as we see from eq. (|A4|) , / n (/?)'s are ordinary numbers and therefore it is not 
possible to satisfy (|A6|). This shows that as long as we restrict ourselves to the original 
Hilbert space, we cannot define a finite temperature vacuum which would satisfy eq. ( |A3|) . 
Therefore we have to double the Hilbert space, introducing a fictitious system identical to 
the original system. Let us denote this auxiliary system as a tilde system. Mathematically 
this means that 
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\n,m) = \n) ® \m). (A7) 
Now we can expand \0(/3)) in terms of these \n,fh) states 

W)> = E/n(/3)|n,n) = £/„(/?) |n> ® |n). (A8) 

n n 

In this case we note that 

(0(/3)|A|0(/3)) = ^/:(/3)/ m (/3)(n,n|AKm> = 

n,m 

= E /:(/3)/m(/5)(n|A|m)5 nm = Y,f:(P)W)(n\A\n). 



(A9) 



Here we have used the fact that an operator of the physical system does not act on the 
Hilbert space of the tilde system and vice versa. So, using the orthonormality of the states 
we can write 



(A10) 



(n, m\A\n', m!) = 


{n\A 


n')(fh 


m') 


&mm' 


{n\A\ 


(n, fh \A \n', fh') = 


(n\A 


n')(fh 


fh!) 


3mm' 


(n\A\ 



^From ( |A3| ) and ( |A6|) it is seen that 

fnW) = W) = Z- 1 ' 2 ^) exp(-(3E n /2) (All) 

Now we can write the Fock space based on the vacuum |0(/3)) which is defined by a 
Bogoliubov transformation as 

W)) = U(0)\O,O) (A12) 

where 

U{(3) = exp[-6(aa - altf)] (A13) 



where cosh(#) = 1/y (1 — exp(— (3uj)) and to is the energy. Then the annihilation operators 
at finite temperature are defined as 

a(/3) = U-\P)aU(P), a(J3) = U~\(3)aU{f3) (A14) 

such that 
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a(/3)|0(/3)) = , a(/?)|0(/?)> = 0. (A15) 
Then Fock space is given as 

|0(/3)), at(/3)|0(/3)), at(/3)|0(/3)),... 



(A16) 

^^[at(/3)]«[St(/5r|0(/?)) 

where n and m are the number of physical particles and tilde particles respectively. This 
therefore shows that we can introduce a temperature dependent vacuum state so that the 
statistical ensemble average of any operator can be identified with the expectation value 
of the operator in this state. This, however entails a doubling of the Hilbert space. The 
advantage, on the other hand, lies in the fact that the techniques of zero temperature field 
theory can now be carried over to the finite temperature case. 



Here are the tilde conjugation rules |M| 

(AB) = AB, 

(c x A + c 2 B) = c\A + c* 2 B, 
(it) = (i)t, 
|0C9)> = |0(/3)). 

For example the tilde partner of the weak interaction Lagrangian 



(A17) 



L w = ^ M r a /ierV e (A18) 



(in usual notations) is 



L w = ^ e f Q g^f%, (A19) 

but not —^=i / fl T a fj J eT a i' e which would give zero contribution to the amplitude of the process 
ji — > ev e v^ at tree level . 



B. Green's functions. 

In TFD the thermal doublet for each field is defined as 
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da) I 



\ m 1 






I m J 




i M%) j 



*) = ; = (A20) 





H 




\ i j 




[ M%) J 



^(x) = { y } = { - ' ) (A21) 



where t means transposition with respect the spinor index, and a(= 1, 2) specifies a compo- 
nent of the thermal doublet. The first component is physical, and the second is fictitious. 
The scalar free field may be written as |J: 



(X) = 77r~TT77T f ^fcwT^ + 



A Fermion field may be written as : 

^(x) = W rfp7V p [ CpisM (p, s )e-^ x +4^(p, S )e^] 

s J 

$(x) = W dpN P H, s u(p, s)e'P x + d p , s v(p, s)e- i P x ] 
^( x ) = J2j dpN p [c\ iS U{p : s)e~ wx +d PiS d(p,s)e 



(A22) 



(A23) 



tpx-\ 



where N p = \jraj (2tt) 3 E p , E p = y m 2 + (p) 2 , £t(p, s) = u*(p, s) etc. 

The boson and fermion propagators, A &(A;) and S a b(k) respectively, are 2 ® 2 matrices. 
They are defined as: 

iA ab (x-y) = (O(/3)|T[0 o (x)0 6 (y)]|O(/3)) =i [ ^e~^ x ~ v)A° b {q) 

7 W (A24) 
iS o6 (* - y) = (0(/3)|T[^(x)^(y)]|0(/3)) = < / -^f~ iq{x ~ ^V<?) 

The free propagators may be presented as a sum of Feynman propagator, A° F (k) (S F (k)) 
and temperature dependent part, A T (A;), (S T (k)) as: 
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A° u (k) = A° F (k) + A§.(fc), A° 2 (k) = -(&°n(k)y 
A° 12 (k) = -2m5{k 2 - m 2 )m B {k), A° 21 (k) = A° 12 (k) 

&° F (k) = j- 2 \ A%{k) = -2mn B (k)5(k 2 - m 2 ), 

k 2 — m 2 + te 

S° u (k) = {Jk + m)(S° F (k) + SUk)), S° 22 (k) = {jk + m)(S° F (k) + S°(k))* (A25) 

S® 2 (k) = 2m{Jk + m)5{k 2 - m 2 )m F (k) S Q 2l (k) = S° l2 {k) 
1 



S F (k) — — 

k 2 — m 2 + ie 

S°(k) = 2iir5{k 2 -m 2 )sm 2 6 h) . 
Here the following notations are used 

sm 2 6 h) = Q(k )n F (k) + Q(-k )n F (k), 

= W 2 [e ( y-e(-M] = f^M sin 20+k ( A26) 

e x + 1 2 
m B (A;) = eP k o/ 2 n B (k) = §sinh20 fe , x = (3{k - fi) 

where 9 is the step function : 6(/c ) = 1 if k > and Q(k ) = otherwise and sign(k ) is 
the sign of ko- 

The full propagator satisfies the Dyson equation: 

^ ab (k) = A° ab (k) + j2 K c m? d (k)A db (k) 

cd (A27) 
Sab(k) = S° ab (k) + £ Sl{k)H F cd {k)S db {k) 

cd 

where Xf b (A:) and are the proper self energies of a fermion and boson respectively. 

By introducing a complex function E B = i?eE B + UmE B we can represent the general form 
of the self energy as follows: 



Ef^fc) = ReE B (k) + UmZ B (k) cosh20 fc , Ef 2 (fc) = -(Ef^/c))" 
Ef 2 ( fc) = -i sinh 26 k Imt B {k) , E^ (Jfe) = Ef 2 (Jfe) . 
Similarly for the fermion self energy we have: 

Ef 1 (/c) = ReTj F (k) + i cos29 k ImTj F (k) H^ 2 (k) = (££(*;))* 
Ef 2 (A;) = -2iImZ F (k)m F (k), E^(ifc) = Ef 2 (ib). 
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(A28) 



(A29) 



The full propagator has the following compact form in terms of E: 

cosh 2 ^ sinh 2 6» fc 



k 2 — m 2 — ReTisik) — HmT, B (k) + is k 2 — m 2 — i?eEs(A;) + UmTtsik) — is 



A 12 (k) m B {k){ k2 _ m2 _ ReTiB ( k ) _ U m E B (k) + is k 2 - m 2 - ReZ B {k) + z/mE B (A;) - ie' 



A 22 {k) = -(A 11 (k))*, A 21 {k) = A 12 (k) 

(A30) 



Sn{k) 



cos 2 6k sin 2 9 ko 



m — ReT l p(k) — ilmT^p^k) + is fa — m — ReHp{k) + HmT,p(k) — is 



^ F ^ ^ fa — m — ReYip{k) — HmY>p{k) + ie fa — m — ReT,p{k) + HmEp(k) — ie^ 



S 21 {k) = S 12 {k) S 22 (k) = (Su(k))* , 

(A31) 

where cosh 2 9 k = 1 + n B , sinh 2 9 k = ns{k), cos 2 9k — 1 — sin 2 and sin 2 9k was defined 
in Eq. ( |A26| ). The spectral representation for propagators and their matrix form may be 
found in ref. 01 . 
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iA+(q) iA-(q) 




iA (p — q) iA + (p — q) 



Fig. 1. Generalized Cutcosky rules for the Imaganary part of 
boson self energy. See ref. [11] for the details. 
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Fig. 2 Feynman diagrams for boson - boson scattering amplitude 
up A 2 in X(j) 4 (x)/4\ model at finite temperature. Here the following 
notations are used : C(k) = cosh^, S(k) = sinh^ , K = k\ + k 2 = 
ki + k 2 ' and A e = X 2 N £ /16n 2 
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